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PREFACE 


This  report  describes  a method  that  greatly  facilitates 
the  analysis  of  many  factorial  designs  and  is  expected  to  be 
especially  helpful  in  developing  a catalogue  of  fractional 
factorial  experiment  designs  for  use  in  experiments  in 
which  some  factors  are  to  be  studied  at  two  conditions  and 
others  at  three  conditions.  This  work  is  under  the 
immediate  direction  of  W.  S.  Connor  and  M.  Zelen.  Professor 
R.  C.  Bose  of  the  University  of  North  Carolina  is  serving  as 


consultant. 
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Introduction , 

The  problem  of  constructing  fractional  factorial 
designs  for  the  2m3n  series  has  not  received  much  attention 
in  the  statistical  literature  (see,  e.g.,  [l]).  Accordingly, 
before  preparing  a catalogue,  it  is  desirable  to  study 
methods  of  construction. 

Such  a study  is  now  in  progress.  The  first  major 
problem  encountered  was  the  formation  of  the  normal  equations 
which  correspond  to  any  particular  design.  Although  it  is 
known  from  the  general  theory  of  least  squares  how  to  form 
the  equations  of  expectation,  and  from  them  the  normal 
equations,  we  were  still  faced  with  the  difficult  problem 
of  implementing  thq  theory.  We  recognized  that  progress 
would  be  very  slow  indeed  unless  some  short  cuts  could  be 
found . 

Fortunately,  we  have  found  a very  effective  method 

for  forming  the  normal  equations.  The  method  is  described 

in  this  report,  and  is  applied  to  a 1/2  replicate  of  the 
3 2 

2 3 design.  The  use  of  this  fractional  replicate  is 
illustrated  by  application  to  some  real  experimental  data. 

The  method  is  so  general  that  it  applies  to  the  analysis 
of  all  factorial  experiments.  It  can  be  used  for  analyzing 
fractional  factorials  of  any  mixed  series  - not  merely  the 
2m3n  series.  Another  use  is  for  complete  or  fractional 
factorials  with  repetitions  of  some  treatment  combinations. 
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We  are  now  engaged  in  using  the  method  to  study  various 
designs  for  several  fractions  of  the  2ra3n  series.  Certain 
considerations  of  symmetry  suggest  that  particular  designs 
are  optimum- in  the  sense  of  providing  estimates  which  are 
not  highly  correlated.  However,  these  estimates  may  be 
arithmetically  more  complicated  than  the  estimates  which  are 
provided  by  other  less  symmetrical  arrangements.  Thus, 
designs  which  are  optimum  in  one  sense  are  not  necessarily 
optimum  in  some  other  sense . 

We  plan  to  program  the  method  for  use  on  the  National 
Bureau  of  Standards  electronic  computer.  We  then  shall  be 
able  to  progress  rapidly  in  the  evaluation  of  the  comparative 
merits  of  various  contending  designs.  From  these  we  shall 
select  "optimum"  designs  for  inclusion  in  the  catalogue. 
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The  method  of  forming  the  design 

We  partition  the  N factors  into  two  collections , and 
C2>  containing  c^  and  c 2 factors,  respectively,  (Ci+C2  “N) , 
Then,  for  each  collection  separately,  we  form  w sets  of 
treatment  combinations.  It  is  not  required  that  all  of  the 
treatment  combinations  be  included  in  some  set,  nor  is  there 
any  restriction  as  to  the  number  of  sets  which  contain  any 
particular  treatment  combination.  Also,  a treatment 
combination  may  occur  more  than  once  in  a set. 

We  denote  the  sets  for  collection  by  m^,  m^,  ..., 
m^  and  their  respective  numbers  of  treatment  combinations  by 


u 


u 


1'  2 9 


u 


w 


Similarly,  we  denote  the  sets  for  collection 


C2  by  n^,  n2>  ...,  n^  and  their  respective  numbers  of 

treatment  combinations  by  v.. , v_,  ...,  v . 

The  design  consists  of  adjoining  every  treatment 

combination  of  set  m^  to  every  treatment  combination  of  set 

n^,  (i-l,  ...,  w) . This  produces  sets  p^  which  contain 

u.v.  treatment  combinations  for  all  N factors, 
l i 

To  illustrate,  let  there  be  N=3  factors  A,  B,  and  a, 

all  having  two  levels « Let  consist  of  A and  B;  and  C2 

consist  of  a.  We  may  choose  the  sets  m.  and  n.  as  follows: 

ii 
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A B 


a 


m 


1* 


0 0 (u  =2)  n 
0 0 1 1 


0 


(vl"1) 


m 


2* 


A B 

0 0 (u  =2)  n 
11  * Z 


a 


0 

1 


(v2=2) 


m3 ' 


A B 

1 0 (u3=1)  n3 


a 


(v3=1) 


The  design  D consists  of  the  following  treatment 


combinations : 


Pl: 


A 

0 

0 


B a 


0 


0 


(ulvl“2) 


P2: 


0 

0 

1 

1 


0 

0 

1 

1 


0 

1 

0 

1 


^U2V2 


V 


0 


(U3V3=1) 
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Equations  of  expectation  and  normal  equations 

We  introduce  column  vectors  of  expected  responses 
(m^)  and  r\(n/) , and  write  the  equations  of  expectation  as 

7)(mi>  = NL 
and  T\(ni)  = Ni 


g 

q 


g 

P 


where  g denotes  the  grand  mean,  p is  a column  vector  of  s 

parameters,  q is  a column  vector  of  t parameters:  and  ML  is 

a u.  X (s+1)  matrix  of  coefficients  and  N.  is  a v.  * (t+1) 
1 11 

matrix  of  coefficients. 

We  now  are  ready  to  consider  the  equations  of 
expectation  and  the  normal  equations  for  the  design  D. 

For  this  purpose,  we  introduce 


y(D) 


y(px) 


_y(pwL 

a column  vector  of  the  responses  of  the  treatment  combina- 
tions of  D,  r^(D) , a column  vector  of  the  corresponding 
expected  responses,  and  r , a column  vector  of  st  parameters. 

We  shall  be  concerned  with  the  following  equations  of 
expectation  for  the  treatment  combinations  of  D: 
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Y|(D)  - 
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M1  a N1 

g 

M2  ® N2 

P 

* 

q 

MSN 

r 

w w 

— 

- __ 
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_ _L 

where  £ denotes  the  right  direct  product*  Hie  normal 
equations  are  as  follows: 


[(M1aN1)T  ,..(MwaNw)T] 


(M10N1) 

— ^ 
g 

• 

P 

q 

(MSN) 

r 

w w 

_ j 

= [(M1SN1)T  ...  (M%aNw)T]  y (D) 


By  examination  of  the  indicated  operations  it  can  be 

verified  that 

(M.SN.)T  = m.tsn.t 
v 1 1 11 


and 


(M.SN.)T  (M.8N.)  = M.TM.8N.TN. 

lx  i,iii 


Hence,  the  normal  equations  become 

T„  T 


[X  (M.  M.SN.  N.)] 

1 i 1 1 1 1 J 


P 

q 

r 


Z(M±®K±)l  y(p.) 


* If  A=*(a.  .)  and  B=(b,  ,),  then  AfiB^a.  .B).  See,  e*g.,[4] 
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De f inition  of  the  parameters . 


The  parameters  g,  p,  q and  r will  be  defined  in  terms 
of  the  expected  responses  in  the  complete  factorial. 

Let  the  equations  of  expectation  for  the  complete 
factorial  (CF)  be  written  as 

r\  (CF)  - E P (CF ) , 

where  (CF)  is  a column  vector  containing  (s+1) (t+1) 
expected  responses,  E is  a square  matrix  of  order  (s+1) (t+1) 
and  p(CF)  is  a column  vector  containing  (s+1) (t+1) 

parameters.  If  E is  non-singular,  then  the  parameters 

\ 

p(CF)  are  defined  by 

£(CF)  - E-1»| (CF) 

We  shall  define  the  matrix  E as  follows.  The  vector 
of  coefficients  of  g (a  column  of  E)  is  the  unit  vector. 

For  any  factor  F with  i levels,  we  introduce  i-1  parameters, 
which  may  be  referred  to  as  the  linear  effect,  the 
quadratic  effect,  the  cubic  effect,  etc.  The  elements 
in  the  vector  of  coefficients  (a  column  of  E)  of  each  of 
these  effects  are  the  values  of  the  corresponding 
orthogonal  polynomial  [2]. 

For  example,  if  a factor  has  three  levels,  0,  1 and  2, 
then  there  are  two  parameters:  the  linear  effect  and  the 

quadratic  effect.  The  coefficient  of  the  linear  effect  is  -1 
0,  or  1 and  of  the  quadratic  effect  is  1,  -2  or  1,  depending 
on  whether  the  factor  is  at  level  0,  1 or  2 . 
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Other  parameters  are  associated  with  several  factors 
simultaneously,  and  may  be  called  interaction  effects.  The 
vector  of  coefficients  of  such  an  effect  (a  column  of  E) 
is  obtained  by  taking  the  element  by  element  products  of  the 
several  corresponding  vectors.  For  example , consider  factors 
Fp  and  F^  which  have  2,  2 and  3 levels,  respectively. 

The  vector  of  coefficients  of  the  linear  by  linear 
by  quadratic  effect  is  obtained  as  indicated  below: 

Level  Effect 

F^  Linear ,F^  Linear ,F^  Quad. ,F^  L*L*Q 


0 

0 

0 

-1 

-1 

1 

1 

0 

1 

0 

-1 

1 

1 

-1 

1 

0 

0 

1 

-1 

1 

-1 

1 

1 

0 

1 

1 

1 

1 

0 

0 

1 

-1 

-1 

-2 

-2 

0 

1 

1 

-1 

1 

-2 

2 

1 

0 

1 

1 

-1 

-2 

2 

1 

1 

1 

1 

1 

-2 

-2 

0 

0 

2 

-1 

-1 

1 

1 

0 

1 

2 

-1 

1 

1 

-1 

1 

0 

2 

1 

-1 

1 

-1 

1 

1 

2 

1 

1 

1 

1 

We 

have  defined  E, 

and  now  shall 

examine  the 

def ini 

of 

the 

parameters . 

We 

may  write 

E (dl)  = C 

where  (dl)  is  a diagonal  matrix  and  C is  an  orthogonal 

2 -1/2 

matrix.  The  element  d.  . of  (dl)  is  (2  e ) , where  e.  . 

33  i ij  13 

is  the  element  in  the  ith  row  and  jth  column  of  E. 
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Because  E(dl)  is  orthogonal, 

[E (dl) ]_1  - [E(dl)f, 

whence  9 T 

E = (dl)  E 

Therefore , 

p (CF)  - (dI)2ET  jg(CF)  . 

This  equation  defines  the  effects  in  terms  of  the 
expected  responses. 
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3 2 

A one-half  replicate  of  the  2 3 . 

The  following  set  of  tables.  Tables  1,  . ..,  8 traces 

the  application  of  the  method  to  a particular  1/2  replicate 
3 2 

of  the  23  complete  factorial. 

Table  1 contains  the  36  treatment  combinations  which 


comprise  the  design,  together  with  data  for  an  illustrative 
example.  There  are  three  factors.  A,  B,  and  C which  have 
two  levels,  and  two  factors,  a and  p,  which  have  three 


levels.  The  collection  contains  A,  B,  C and  the  collection 


Cg  contains  a and  p.  The  sets  of  are 


ABC 
0 0 0 
110 
10  1 
Oil 


and  m^ : 


ABC 
10  0 
0 10 
0 0 1 
111 


and  the  sets  of  C^  are 

a J3 
0 0 
0 2 

n^:  11  and  n^ : 

2 0 
2 2 


a p 
0 1 
1 0 
1 2 
2 1 


The  treatment  combinations  of  m^  are  adjoined  to  those 
of  n, , and  those  of  m to  those  of  n to  form  the 

X 4 ^ 

4*5  + 4*4  = 36  treatment  combinations  in  Table  1. 
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TABLE  1 


Treatment  Combinations* 
(Observed  values  in  parentheses) 


A 

0 

1 

1 

0 

1 

0 

0 1 

0 

110 

B 

0 

1 

0 

1 

0 

1 

0 1 

0 

10  1 

C 

0 

0 

1 

1 

0 

0 

1 1 

0 

Oil 

a 

0 

0 

0 

0 

0 

0 

0 0 

0 

0 0 0 

P 

0 

0 

0 

0 

1 

1 

1 1 

2 

2 2 2 

Z"N 

ZN  ZN  z"V 

/*-N 

05 

/•s 

lO 

Z— v 

oo 

z~\ 

CO 

z~v 

o 

/ — > 
rH 

z-N  Z-N 

05  o 

00 

• 

CO  05  CM 
# • • 

lO 

00 

in 

rH 

rH 

05 

rH 

rH 

05 

05 

• 

CM 

rH 

• 

oo 

• • 
00  CM 
00  rH 
rH 

rH 

CO 

rH 

O CO  CO 
05  o CM 
H CM  rH 
sz  v sz 

s»z 

N-Z 

v-/ 

V-z 

s-z 

A 

1 

0 

0 

1 

0 

1 

1 0 

1 

0 0 1 

B 

0 

1 

0 

1 

0 

1 

0 1 

0 

10  1 

C 

0 

0 

1 

1 

0 

0 

1 1 

0 

Oil 

a 

1 

1 

1 

1 

1 

1 

1 1 

1 

111 

P 

0 

0 

0 

0 

1 

1 

1 1 

2 

2 2 2 

Z"N 

z-\ 

z— s 

z-N 

z— \ 

z— s 

z*N  z-N 

z-N 

Z-N  z-\ 

00 

rH 

05 

CM 

00 

CM 

CO  rH 

CM 

CO  rH  CO 

CM  O 

rH 

l> 

rH 

00 

rH  rH 

O 

05  rH  in 

05  iH 

05 

CO 

m 

CM 

rH  rH 

rH 

05  o o 

rH 

rH 

rH 

rH  rH 

rH 

rH  CM 

'w' 

V-z 

w 

w 

v-z  W 

VZ  V— z vz 

A 

0 

1 

1 

0 

1 

0 

0 1 

0 

110 

B 

0 

1 

0 

1 

0 

1 

0 1 

0 

10  1 

C 

0 

0 

1 

1 

0 

0 

1 1 

0 

Oil 

a 

2 

2 

2 

2 

2 

2 

2 2 

2 

2 2 2 

P 

0 

0 

0 

0 

1 

1 

1 1 

2 

2 2 2 

z-s. 

z-s 

Z—N 

/■N 

Z-N  ^ 

Z"N 

Z-N  ZN  Z*S 

CO  CM 

m 

05 

co 

O CM 

O in  rH 

© O 

oo 

m 

rH 

t> 

rH  05 

00 

CO  CM 

CO 

rH 

rH 

CM 

CM 

rH  00 

CO  05  t> 

V-Z 

rH 

S-Z 

rH 

rH 

rH 

w 

rH 

v-z 

rH  rH 

N-Z  V-Z 

rH 

v-z 

rH  rH  rH 
s-'  ^ v— z 

* The  entries  are  read  vertically,  in  groups  of  five. 
For  example,  the  first  entry  indicates  that  all  of  the 
factors  are  at  their  zero  level. 
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TABLE  2 

Expected  Values  of  the  Responses 
iti  a 23  Factorial 


Effect 


Response  at 
Combination 


000 

110 

101 

Oil 


o 


+ 


O 

o 

* 

O 

o 

/*N 

O 

o 

<0 

o 

H 

O 

s 

o 

H 

H 

H 

3 

f? 

S 

M 

(4 

M 

H 

- 

+ 

+ 

4* 

4* 

PM 

+ - + 


+ 4* 


100 

010 

001 

111 


+ 

m 

+ 


4- 


4- 


4* 


+ 

4» 


B(111,00) 
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table  3 

Expected  Values  of  the  Responses 

o 

in  a 3 Factorial 


Effect 


Response  at 
Combination 


00 

02 

11 

20 

22 


01 

10 

12 

21 


to  HI 


24 


TABUS  4 

The  Matrices  M^M^  and 


mJmi 


0 

4 


0 

0 

4 


0 

0 

0 

4 


Symmetric 


0 

0 

0 

-4 

4 


0 

0 

*•4 

0 

0 

4 


0 

-4 

0 

0 

0 

0 

4 


.4 

0 

0 

0 

0 

0 

0 

4 


“2 


0 

4 


0 

0 

4 


0 

0 

0 

4 


Symmetric 


0 

0 

0 

4 

4 


0 

0 

4 

0 

0 

4 


0 

4 

0 

0 

0 

0 

4 


4 

0 

0 

0 

0 

0 

0 

4 


15 


»i«i 


NI»2 


TABLE  5 

The  Matrices  and  m|»2 


5 0 2 0 

4 0 0 

8 0 
4 


Symmetric 


2 0 0 

0 0 4 

0 0 0 

0 0 0 

8 0 0 

4 0 

4 


4 0 -2  0 -2 

2 0 0 0 
10  0 -8 

2 0 
10 


Symmetric 


0 0 

0 -4 

0 0 

0 0 

0 0 

0 0 

8 


0 8 
0 0 
0 -4 

4 0 

0 -4 

0 0 
0 0 
4 0 

20 


0 -8 
0 0 
0 4 

-4  0 

0 4 

0 0 
0 0 
S 0 
16 


TABLE  6 

The  Coefficients  in  the  Normal  Equations 


(mJmj  e nJm^)  + (mJm2  e hJ»2) 


<?X<7 


17- 


and  B * 


-4 (W^NX 


*J»2) 


0 

2 


4 

0 

•2 


0 

0 

0 

2 


Symmetric 


4 

0 

19 

0 

-2 


0 

6 

6 

b 

b 

4 


0 

8 

0 

0 

0 

0 

-4 


0 

0 

0 

8 

0 

0 

0 

•-4 


16 

0 

-8 

0 

-8 

0 

0 

0 

4 


c){ci 
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TABLE  7 

Effects  Which  Occur  in  Each  Subset  of  the  Normal  Equations 


Subset 


1 

2 

3 

4 

E (000, 00)* 

B (100, 00)* 

E (010, 00)* 

E (001, 00)* 

E (000, 10)* 

E(100, 10)* 

E (010, 10)* 

£(001,10)* 

E (000, 20)* 

E(100, 20)* 

E (010, 20)* 

E (001, 20)* 

E(000, 01)* 

E (100, 01)* 

E(010,01)* 

E (001, 01)* 

B (000, 02)* 

E (100, 02)* 

E (010, 02)* 

E (001, 02)* 

E (000, 11)* 

E(100, 11) 

E(010, 11) 

E (001, 11) 

B(000, 12)* 

8,(100, 12) 

E (010, 12) 

E (001, 12) 

E (000, 21)* 

E (100, 21) 

E(010,21) 

E (001, 21) 

B(000,22)* 

B (100, 22) 

E (010, 22) 

E(001,22) 

E (111, 00) 

E(011,00)* 

E(101, 00)* 

£(110,00)* 

B(lll, 10) 

E (011, 10) 

E (101, 10) 

E(110, 10) 

E (111,20) 

E(011,20) 

E(101,20) 

£(110,20) 

E(lll,01) 

E (Oil, 01) 

E (101, 01) 

£ (110, 01) 

E(lll, 02) 

E (Oil, 02) 

E (101, 02) 

E (110, 02) 

E (111, 11) 

8(011,11) 

E(101, 11) 

8(110,  ID 

E(lll, 12) 

E(011, 12) 

E (101, 12) 

£ (110, 12) 

E(lll, 21) 

E(011, 21) 

E(101,21) 

£(110,21) 

E (111, 22) 

8(011,22) 

E (101, 22) 

E(110,22) 

*)  Indicates  the  grand  mean,  the  main  effects,  and  the  two 
factor  interaction  effects* 
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table  a 


Estimation  ol  the  Main  .and  Two -factor  Interaction 

Effects  in  Subset  2 


9 -1 

1(100,00) 

Y (100, 00) 

6 ° 

E (100, 10) 

Y (100, 10) 

18  v_y  -4 

3(100,20) 

tm 

Y (100, 20) 

C)  6 ° 

E (100, 01) 

Y (100, 01) 

V ' 18  -4 

E(100,02) 

Y (100, 02) 

_ -1  0 -4  0 -4  9_ 

_E  (Oil,  002. 

Y (011, 00) 

1 

'7? 


"1(100,00") 

"65  0 2 0 

2 

+{T 

“Y  (100,  00)” 

E (100,10) 

96  0 0 

0 

0 

Y (100, 10) 

E(100, 20) 

* 1 

36  0 

4 

418 

Y (100, 20) 

3 (100,01) 

2 B32 

96 

0 

0 

Y (100, 01) 

E (100, 02) 

36 

+18 

Y (100, 02) 

_1 (Oil, 002 

Symmetric 

81 

Y (Oil, 00) 

, . til  i ■■  ■■■■  •■■■■-• 


' 

' 
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Tables  2 and  3 contain  sketchs  of  the  equations  of 

expectation.  : At  the  left  are  the  expected  responses 

^(000),  r[  (110) , etc.,  with  Y[  ( ) omitted.  The  equality 

signs  have  been  omitted  and  the  parameters  (effects)  have 

been  indicated  at  the  top,  instead  of  as  column  vectors. 

In  Table  2 it  is  understood  that  the  + and  - are  to  be 

read  as  +1  and  -1.  The  matrices  M., , M . N,  and  have 

1’  2 1 2 

been  indicated. 

The  notation  for  the  effects  is  as  follows:  The  symbol  E 

denotes  "effect”,  and  the  five  positions  correspond,  in  order, 
to  A,  B,  C,  a and  (3.  The  comma  separates  A,  B,  and  C from 
a and  (3.  Zero  (0)  in  a position  indicates  that  no  effect 
of  the  corresponding  factor  is  involved  in  the  definition 
of  S(  ).  One  (1)  in  a position  indicates  that  the 

linear  effect  of  the  corresponding  factor  is  involved,  and 
two  (2)  that  the  quadratic  effect  is  involved. 

The  following  typical  examples  should  suffice  to 
explain  the  notation: 


Grand  average 
Main  effect  of  A 
Interaction,  A with  B 
Linear  effect  of  a 
Quadratic  effect  of  a 
Linear  a by  linear  (3 
Linear  a by  quadratic  (3 
Linear  A by  linear  a 
Linear  A by  quadratic  a 


E (000,00) 
E(100,002 
E (110,00) 
E (000, 10) 
E (000,20) 
E (000, 11) 
E (000, 12) 
E (100, 10) 
E (100,20) 


T T 

Tables  4 and  5 exhibit  the  matrices  M^,  M2  M^, 

T T 

N..  N, , and  N NQ  ; and  Table  6 displays  the  matrix  of 
1 i £ £ 

coefficients  of  the  normal  equations.  This  matrix  is  of 
order  72 , and  decomposes  into  four  identical  submatrices 
of  order  18.  These  submatrices  are  of  rank  9,  which  implies 
that  it  is  impossible  to  uniquely  estimate  the  72  effects. 

At  this  point  it  is  reasonable  to  equate  certain  effects 
to  zero,  and  to  solve  for  the  remaining  effects.  A selection 
of  effects  to  be  estimated  is  indicated  by  asterisks  in 
Table  7.  (The  unstarred  effects  are  put  equal  to  zero.) 

The  effects  to  be  estimated  are  all  of  the  main  effects 
and  two-factor  interactions . 

The  subsets  in  Table  7 correspond  to  the  submatrices 
of  Table  6.  Solutions  for  the  effects  in  subset  1 are 
particularly  easy,  because  the  corresponding  matrix  of 
coefficients  is  a diagonal  matrix. 

Table  8 displays  the  solution  for  the  effects  in 
subsets  2,  3 and  4.  The  symbol  YCi^igi^,  ^1^2^  ~ where 
il,  i2>  and  i^  take  on  the  values  0 and  1;  and  and 
take  on  the  values  0,  1 and  2 - denotes  a linear  combination 
of  the  observed  responses.  It  occurs  in  the  normal  equation 
which  corresponds  to  ECi^^ig,  and  the  ianer 

product  of  the  vector  of  responses  with  the  vector  of 
coefficients  of  E(i  i i , j-,j0)  in  the  equations  of 

I Cl  «3  I 


expectation . 


From  the  inverse  matrix,  one  may  read  the  variances  and 
covariances  of  the  estimates.  The  maximum  correlation 
between  estimates  is  18/[  (36)  (81)  = 1/3. 

There  is  very  little  correlation  among  the  estimates. 

In  fact,  of  the  27(26)/2  = 351  pairs  of  estimates,  there  is 
correlation  between  the  estimates  in  only  18  pairs. 

If,  instead  of  using  the  present  design,  one  used  the 
design*which  is  formed  by  adjoining  the  treatments  of  m^  to 
those  of  n,  and  n , only  21  of  the  chosen  effects  would  be 

* Z 

estimable  (the  effects  E (100,00),  E (Oil, 00),  E (010, 00), 

E (101, 00),  E (001, 00),  and  E (110,00)  would  not  be  estimable). 

It  is  interesting  to  compare  the  variances  of  the 
estimates  obtained  from  the  two  designs.  They  are  identical 
except  in  six  instances.  For  E (100, 20),  E (100, 02), 

E (010,20) , E (010,02) , E(001,20),  and  E(001,02)  the 
variance  for  the  alternative  design  is  8/9  as  large  as 
for  the  present  design. 


* As  suggested  in  [1]. 
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An  illustrative  example 


The  following  example  is  taken  from  a paper  by  Youden  and 
Zimmerman  [3].  Youden  was  concerned  with  comparing  various 
methods  of  producing  tomato  plant  seedlings  prior  to  trans- 
planting in  the  field.  Comparison  was  made  by  planting  in 
the  field  and  then  weighing  the  ripe  produce.  Thus,  the 
observations  were  pounds  of  tomatoes. 

Although  Youden  used  five  methods  of  production,  we  shall 
select  only  three:  flats,  fibre  pots,  and  fibre  pots  soaked 

in  one  percent  sodium  nitrate  solution.  Other  factors 
considered  were  different  soil  conditions,  different  sizes 
of  pots,  different  varieties  of  tomato,  and  different 
locations  on  the  field.  The  factors  and  their  levels  are 
recorded  below: 


Factor 

Soil  condition,  A 
Size  of  pot,  B 
Variety  of  tomato,  C 
Method  of  production,  a 


Levels 

Field  Soil  0 

Plus  fertilizer  1 

Three-inch  0 

Four- inch  1 

Bonny  Best  0 

Margiobe  1 

Flat  0 

Fibre  1 

Fibre  + N0„  2 


0,1,2 


Location  on  field,  £ 
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The  object  of  the  experiment  was  to  evaluate  the 
effects  of  these  factors  on  the  yield  of  tomatoes.  The 
observations  have  been  recorded,  along  with  the  experiment 
design,  in  Table  1. 

The  first  stage  in  the  analysis  is  the  same  as  for  a 
complete  factorial.  We  form  a number  of  summary  tables, 
which  are  given  below: 

TABLE  9 

Summary  Tables 


Soil 

Condition 

(A) 

0 

1 

Total 

Size  of 

0 

993.7 

1210.9 

2204.6 

Pot  (B) 

1 

1100.7 

1415.9 

2516.6 

Total 

2094.4 

2626.8 

4721.2 

Soil 

Condition 

(A) 

0 

1 

Total 

Variety 

0 

980.3 

1078.8 

2059.1 

(C) 

1 

1114.1 

1548.0 

2662.1 

Total 

2094.4 

2626.8 

4721.2 

Size 

of  Pot  (B) 

0 

1 

Total 

Variety 

0 

931.1 

1128.0 

2059.1 

(C) 

1 

1273.5 

1388.6 

2662.1 

Total 

2204.6 

2516.6 

4721.2 

Method  of  Production  (a) 

0 

1 

2 

Total 

Soil  Cond. 

0 

640.5 

608.2 

845.7 

2094.4 

(A) 

1 

813.5 

848.3 

965.0 

2626.8 

Total 

1454.0 

1456.5 

1810.7 

4721.2 
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Table  9 (Continued) 

Method  of  Production  (a) 


0 

1 

2 

Total 

Size  of 

0 

705.3 

601.7 

897.6 

2204.6 

Pot  (B) 

1 

748.7 

854.8 

913.1 

2516.6 

Total 

1454.0 

1456.5 

1810.7 

4721.2 

Method  of  Production  (a) 

0 

1 

2 

Total 

Variety 

0 

676.9 

595.7 

786.5 

2059.1 

(C) 

1 

777.1 

860.8 

1024.2 

2662.1 

Total 

1454.0 

1456.5 

1810.7 

4721.2 

Location  (p) 

0 

1 

2 

Total 

Soil  Cond. 

0 

617.0 

635.2 

842.2 

2094.4 

(A) 

1 

784.0 

767.6 

1075.2 

2626.8 

Total 

1401.0 

1402.8 

1917.4 

4721.2 

Location  (p) 

0 

1 

2 

Total 

Size  of 

0 

652.5 

592.9 

959.2 

2204.6 

Pot  (B) 

1 

748.5 

809.9 

958.2 

2516.6 

Total 

1401.0 

1402.8 

1917.4 

4721.2 

Location  (p) 

0 

1 

2 

Total 

Variety 

0 

595 . 4 

552.7 

911.0 

2059.1 

(C) 

1 

805.6 

850.1 

1006.4 

2662.1 

Total 

1401.0 

1402.8 

1917.4 

4721.2 

Method  of  Production  (a) 

0 

1 

2 

Total 

Location 

0 

420.5 

465.3 

515.2 

1401.0 

O) 

1 

351.3 

472.0 

579.5 

1402.8 

2 

682.2 

519.2 

716.0 

1917.4 

Total 

1454.0 

1456.5 

1810.7 

4721  2 
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The  entry  in  any  cell  is  the  sura  of  the  observations  which 
have  the  row  and  column  designations  associated  with  the  cell. 
For  example,  in  the  last  table,  the  first  entry  is  the  sum 
of  the  observations  in  Table  1 which  have  a - 0 and  p * 0, 
i.e.,  the  first  four  observations. 

We  calculate  the  Y's  from  the  summary  tables.  For 
example,  from  the  first  summary  table  we  compute 

Y (100,00)  - -993.7  - 1100.7  + 1210.9  + 1415.9  =■  532.4 

Y (010,00)  993.7  + 1100.7  - 1210.9  + 1415.9  = 312.0 

and  Y (110,00)  = 993.7  - 1100.7  - 1210.9  + 1415.9  = 98.0 

We  observe  that  there  is  some  repetition  of  Y's  from  the 

various  tables.  For  instance,  Y (100, 00)  is  obtained  from  all 

four  tables  which  have  A as  a factor.  This  repetition  is 

unnecessary  and  in  practice  would  be  avoided.  Having  already 

calculated  a Y from  a previous  table,  we  would  not  calculate 

it  again. 

The  complete  list  of  27  distinct  Y*s  is  given  in  Table  10 


TABLE  10 
Values  of  Y's 


Y1 

= 

Y (000,00) 

= 

4721.2 

Y15 

— 

Y (001,20) 

= 

-192.3 

Y2 

= 

Y (100,00) 

* 

532.4 

Y16 

= 

Y (000,01) 

= 

516.4 

Y 

3 

= 

Y (010, 00) 

= 

312.0 

Y17 

= 

Y (000,02) 

3 

512.8 

Y4 

= 

Y (110,00) 

= 

98.0 

Y18 

= 

Y (100,01) 

3 

66.0 

Y5 

3 

Y (001,00) 

= 

603.0 

Y19 

= 

Y(100,02) 

— 

135.2 

Y6 

= 

Y (101,00) 

= 

335.4 

Y20 

= 

Y (010,01) 

3 

-97.0 

Y7 

= 

Y (011,00) 

as 

-81.8 

Y21 

5 

Y (010,02) 

3 

-339.0 

Y8 

= 

Y (000, 10) 

= 

356.7 

Y22 

= 

Y (001,01) 

3 

-114.8 

Y9 

=1 

Y (000,20) 

- 

351.7 

Y23 

- 

Y (001,02) 

- 

-289.2 

Y10 

3 

Y (100,10) 

3 

-53.7 

Y24 

Y (000, 11) 

-60.9 

Y11 

= 

Y(100, 20) 

- 

-187.9 

Y25 

- 

Y (000,21) 

- 

354.7 

Y12 

- 

Y (010, 10) 

- 

-27.9 

Y26 

- 

Y(000, 12) 

- 

-327.9 

Y13 

3 

Y(010,20) 

- 

-447.3 

Y27 

Y (000, 22) 

391.3 

Y14 

= 

Y (001, 10) 

- 

137.5 

The  formulae  and  estimates  of  the  main  and  interaction 
effects  are  given  in  Table  11. 
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TABLE  II 
Estimated  Effects 


E (000,00) 

= 

Y1/36 

= 131.1 

E (100, 10)  - Y1q/12  - 

-4.5 

E (000, 10) 

= 

Yg/12 

= 29.7 

E (100,01)  - Y18/12  - 

5.5 

E (000,01) 

- 

Y16/12 

■ 43.0 

E (010 , 10)  = Y12/12  - 

-2.3 

E (000, 20) 

Yg/24 

- 14.7 

E (010, 01)  - Y20/12  = 

-$.l 

E (000,02) 

= 

Y17/24 

- 21.4 

E (001, 10)  - Y,  ./12  - 

11.5 

E (000,11) 

= 

Y24/8 

=*  -7.6 

E (001, 01)  - Y22/12  - 

-9.6 

E (000,21) 

= 

Y25/16 

- 22.2 

E (000,12) 

Y26/16 

- -20.5 

E (000,22) 

= 

Y27/32 

- 12.2 

E (100,00) 

= 

(65^ 

+ 2Y  + 2Y 

11 

19 

+ 9Y  )/1152  =29.3 

E(100,20) 

- 

(Y2  + 

18Y11  + 2Y19 

+ 9Y?)/  384  = -8.6 

E (100,02) 

= 

(Y2  + 

2Y±1  + 18Y 

19 

+ 9Y?)/  384  = 4.8 

E (011,00) 

- 

(Y2  + 

2Yn+  2Y 

19 

+ 9Y  )/  128 2.4 

E (010,00) 

= 

(65^ 

+ 2Y13  + 2Y 

21 

+ 9Y_)/1152  =18.9 
6 

E (010,20) 

= 

(Y3  + 

18Y  3 + 2Y 

21 

+ 9Y  )/384  =-14.1 

6 

E (010,02) 

= 

(Y3  + 

2Y13  + 18Y21  4 

• 9Y.)/384  9.5 

6 

E (101, 00) 

= 

(Y3  + 

2Y13  + 2Y21 

+ 

9Y.)/128  = +13.7 

6 

E (001,00) 

(65y5 

+ 2Y1C  + 2Y 
15 

23 

+ 9YJ/1152  = 34.0 
4 

E (001,20) 

=3 

(Y5  + 

18Y15  + 2Y23  4 

• 9Y  )/384  = -6.7 

E (001,02) 

= 

(Y5  + 

2Y15  + 18Y23  4 

■ 9Y  )/384  = -10.7 

E (110,00) 

= 

(Y5  + 

2Y15  4 2Y23 

+ 

9Y.)/128  = 4.1 

4 

4.1 
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These  formulae  have  been  arranged  into  four  sets. 
All  of  the  estimated  effects  in  the  first  set  are 


uncorrelated . 

In  addition 

the  estimated 

effects  within 

any  set  are  uncorrelated  with  the  estimated  effects  in 

any  other  set. 

The  analysis  of  variance  is  given  in 

Table  12 . 

Table 

12 

Analysis  of 

Variance 

Source  of 

Degrees  of 

Sum  of 

Mean 

Variation 

Freedom 

Squares 

Square 

Effects 

26 

58975.51 

2268 

Residual 

9 

4150.08 

461 

Total 

35 

63125.59 

In  order  to  evaluate  the  effects,  t-tests  were 
carried  out.  It  was  found  that  E (100, 00),  E (001, 00), 
E(000,10)  and  E(000,01)  are  significant  at  the  .01 
level  of  significance.  There  were  no  other  significant 


effects . 
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Chemistry.  Organic  Coatings.  Surface  Chemistry.  Organic  Chemistry.  Analytical  Chemis- 
try. Inorganic  Chemistry.  Electrodeposition.  Molecular  Structure  and  Properties  of  Gases. 
Physical  Chemistry.  Thermochemistry.  Spectrochemistry.  Pure  Substances. 

Mechanics.  Sound.  Mechanical  Instruments.  Fluid  Mechanics.  Engineering  Mechanics. 
Mass  and  Scale.  Capacity,  Density,  and  Fluid  Meters.  Combustion  Controls. 

Organie  and  Fibrous  Materials.  Rubber.  Textiles.  Paper.  Leather.  Testing  and  Specifica- 
tions. Polymer  Structure.  Plastics.  Dental  Research. 

Metallurgy.  Thermal  Metallurgy.  Chemical  Metallurgy.  Mechanical  Metallurgy.  Corrosion. 
Metal  Physics. 

Mineral  Products.  Engineering  Ceramics.  Class.  Refractories.  Enameled  Metals.  Concret- 
ing Materials.  Constitution  and  Microstruclure. 

Building  Technology.  Structural  Engineering.  Fire  Protection.  Air  Conditioning,  Heating, 
and  Refrigeration.  Floor,  Roof,  and  Wall  Coverings.  Codes  and  Safety  Standards.  Heat 
Transfer. 

Applied  Mathematics.  Numerical  Analysis.  Computation.  Statistical  Engineering.  Mathe- 
matical Physics. 

Data  Processing  Systems.  SEAC  Engineering  Group.  Components  and  Techniques.  Digital 
Circuitry.  Digital  Systems.  Analog  Systems.  Application  Engineering. 

• Office  of  Basic  Instrumentation.  • Office  of  Weights  and  Measures. 

BOULDER,  COLORADO 

Cryogenic  Engineering.  Cryogenic  Equipment.  Cryogenic  Processes.  Properties  ol  Materials. 
Gas  Liquefaction. 

Radio  Propagation  Physics.  Upper  Atmosphere  Research.  Ionospheric  Research.  Regular 
Propagation  Services.  Sun-Earth  Relationships.  VHF  Research. 

Radio  Propagation  Engineering.  Data  Reduction  Instrumentation.  Modulation  Systems. 
Navigation  Systems.  Radio  Noise.  Tropospheric  Measurements.  Tropospheric  Analysis. 
Radio  Systems  Application  Engineering.  Radio  Meteorolog). 

Radio  Standards.  High  Frequency  Electrical  Standards.  Radio  Broadcast  Service.  High 
Frequency  Impedance  Standards.  Calibration  Center.  Microwave  Physics.  Microwave  Circuit 
Standards. 


